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This work presents a step furthering a new perspective of proactive control of the spin-exciton 
dynamics in the quantum limit. Laser manipulation of spin-polarized optical excitations in a semi- 
conductor nanodot is used to control the spin dynamics of two interacting excitons. Shaping of 
femtosecond laser pulses keeps the quantum operation within the decoherence time. Computa- 
tion of the fidelity of the operations and application to the complete solution of a basic quantum 
computing algorithm demonstrate in theory the feasibility of quantum control. 



Experimental laser control of spin states of excitons has 
been demonstrated in ensembles of semiconductor quan- 
tum dots (QDs) 0. In a single dot, ultrafast control of 
spin-excitons , and entanglement of the electron-hole 
complex Q have been reported. Here, we present a study 
of the design of laser pulses for the ultrafast control of the 
spin dynamics of individual excitons in a QD. This is a 
special case of designing Hamiltonians to bring a system 
from one state to another We hope that the theory 
of control of exciton spin dynamics may help rapid re- 
alization of basic quantum operations in nanodots. An 
exciting application would be the further development of 
quantum computation p] using excitons frffol or inter- 




subband transitions pdl-llJ 

In this paper, the control of exciton spin dynamics con- 
sists in designing laser pulses to create as a function of 
time desired multiexciton states of a subsystem of dots. 
In the ultrafast control of the spin dynamics of excitons, 
there is the possibility of unintended dynamics due to the 
presence of close resonances. Sharp resonant pulses can 
minimize such unintended dynamics, but at the cost of a 
long operation time. The latter leads to the spontaneous 
recombination of the excitons and to their dephasing, re- 
sulting in uncontrolled deterioration of the amplitude and 
phase of the coefficients of the desired linear combination 
of quantum states. We give a solution to these two con- 
tradictory requirements by suggesting extending the well- 
known laser pulse-shaping technique Q to the delicate 
control in the quantum limit. Design of excitations has 
been applied to quantum computation by NMR and 
by electron spin dynamics or intersubband transitions in 
QDs We rely on the theory of DiVincenzo ||2|] 

stating that quantum gates operating on just two qubits 
at a time are sufficient to construct a general quantum 
circuit in a system of many dots. As an application to 
a prototype quantum computation which illustrates the 
issues raised above, we present the results of numerical 
simulations of the two exciton dynamics in a sequence 
of operations to solve the Deutsch-Jozsa (DJ) problem 
pT|,p[. We shall concentrate on the control of the lowest 
four states formed by two excitons with opposite spins 



in a single dot. The spin-down and up states can be ex- 
cited by the left (cr— ) and right-handed ((T+) circularly 
polarized light, respectively. The four basis states are, in 
order: the lowest biexciton state | H — ), the a-\- exciton 
state |-|-), the a— exciton state |— ), and the ground state 
0). We use the interaction representation for the states. 
We allow for the presence of other multiexciton states in 
the dot and the connection of the dot to the environment. 

For precise control of coherent exciton dynamics, the 
probabilistic photo-excitation process of an exciton needs 
to be replaced by a resonant Rabi rotation which can cre- 
ate a linear combination of two states at a definite time. 
The primary condition for carrying out a Rabi rotation 
is sufficient light intensity without generating more exci- 
tons. This can be ascertained by analyzing the coherent 
nonlinear optical spectra. The starting point of our cal- 
culation is the energies and wave functions of the multi- 
exciton states in a dot. For the current purpose, it is suf- 
ficient to obtain these states from two confined levels of 
electrons and holes each in a parallelepiped QD . The 
size of the dot, 35 nmx40 nmx5 nm, is typical of inter- 
face fluctuation quantum dots . Only Coulomb inter- 
action between the carriers which conserves their conduc- 
tion or valence band indices is taken into account exactly. 
The resulting states are used to compute the nonlinear 
spectra in a coherent pump and probe experiment setup 
as in |p3|| . Each cw light beam has a narrow frequency 
peak and a definite circular polarization. The two beams 
have a fixed phase difference. The nonlinear spectrum is 
the absorption of the probe beam as a function of its fre- 
quency for a given intensity and frequency of the pump 
beam. An example of the calculated spectrum is shown in 
Fig. |l| for the cross-polarized circular configuration. The 
peaks come from transitions between the multi-exciton 
states. The key to our quantum operations is the lowest 
two peaks, corresponding to the exciton-biexciton and 
the ground state-exciton transitions, respectively (mag- 
nified in the inset). The separation of 1 meV between the 
two peaks corresponds to the binding energy of | H — ). 
An exciton lifetime of 40 ps is taken from the experiment 
p3[ . The doublet structure in each peak is the Rabi split- 
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ting due to the presence of the strong pump field. The 
observation in the nonlinear spectrum of these doublets 
would indicate the feasibility of Rabi rotations in real- 
time operations. The life-time broadening and the line 
positions yield the time limit on the pulses to be used 
in the exciton operations. The observation of Rabi rota- 
tions in QDs has not yet been reported. Evidence of Rabi 
flopping for cxcitons in semiconductor quantum wells p^ ] 
is encouraging. 
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FIG. 1. Coherent nonlinear spectra of a QD with pump and 
probe in opposite circular polarization. The pump is resonant 
with the lowest excitonic state. The solid line is for the pump 
beam with intensity corresponding to a Rabi energy of 50 /^eV 
and the dashed line for zero pump intensity. The inset is a 
zoom on the first two levels. 

The first quantum operation considered here is the con- 
ditional dynamics where the a— exciton will be given a 
TT rotation only if there is already a a+ exciton present. 
This controlled Rabi rotation (C-ROT) is related to the 
controUed-not (C-NOT) or exclusive-or logic operation 
between two qubits Q . A C-NOT operation inverts the 
second qubit (the target) depending on the state of the 
first qubit (the control) . The C-ROT operation is accom- 
plished simply by using a a— circularly polarized laser 
pulse to effect a single resonant Rabi 7r-rotation p5| , p] ,p[ 
which exchanges the exciton state |-|-) with the biexciton 
state — I + — ) [see Fig. ||(a)]. Ideally, the same pulse does 
not cause the transition between |0) and |— ) because the 
interaction between two excitons cause this transition to 
be off-resonance by an amount equal to the binding en- 
ergy of the biexciton. The second operation is a single 
qubit operation: it acts on an exciton regardless of the 
state of the other exciton. It requires a resonant Rabi 
rotation for each of the many transitions involving the 
exciton. In the 2-exciton subspace, an operation on, say, 
the a— exciton requires two frequencies for the two tran- 
sitions between the ground state jO) and |— ) and between 
|-|-) and I -I- — ), i.e., a two-color laser pulse is needed. The 



idea of a resonant Rabi rotation for C-NOT was proposed 
by Lloyd and by Barenco et al, . To complete the 
C-NOT, we could add to a C-ROT additional operations. 
However, we have found no need in quantum algorithms 
to take this time consuming extra step. Proper rotations 
represent the natural choice for the physical design of al- 
gorithms using Rabi flopping. We will show below how 
C-NOT and Hadamard gates in the standard implemen- 
tation of the D J algorithm ||^ can be replaced by C-ROT 
and 7r/2 operations. Using proper rotations, we have also 
designed the optical control for quantum Fourier trans- 
form, which is key to several algorithms including factor- 
ization, and carried out the simulation ||2^ . 

We have simulated these two basic operations and 
calculated the gate fidelity defined in Ref. (2^, whose 
equally weighted average over all the possible initial 
states gives E^_l^«l' + ^ E^^Ji^^IJ3 + = 
{i\U''U\j), where U and U are the calculated and the ideal 
quantum transformations respectively. We have included 
in addition to the four states in the computational basis 
the two copolarized biexciton states | -I- | ) in or- 
der to check their contribution to the unintended dynam- 
ics. The other states are more distant in energy and can 
be neglected. The dephasing is taken into account phe- 
nomenologically as a stochastic process using the Monte 
Carlo wavefunction method |2^ . A 40 ps exciton lifetime 
is taken from j23j. If we use a spectrally sharp pulse to 
perform C-ROT, without dephasing the unintended per- 
turbation is negligible but the operation time is of the 
same order as the dephasing time. On the other hand, 
a sufficiently short pulse (about 2 ps) to limit dephasing 
destroys the fidelity by unintended dynamics. We pro- 
pose a remedy by pulse shaping. Let us first compare 
the results in fidelity of the operations before present- 
ing the details of shaping. The fidelity decreases with 
increasing Rabi frequency. The calculated results used 
2 meV Rabi energy which is larger than the separation 
between the exciton and biexciton resonances of 1 meV. 
The fidelity for C-ROT is 0.535, 0.966, and 0.903 for, in 
order, unshaped short pulse (without dephasing), shaped 
pulse without and with dephasing. For single qubit 7r/2- 
rotation, it is 0.678, 0.995 and 0.981 respectively. The 
fidelity is nicely restored by pulse shaping. The leakage 
to the unwanted states] -I- | ) is less than 1%. 

In a C-ROT, instead of a spectrally narrow tt pulse 
resonant with the exciton-biexciton transition, we use a 
combination of two phase-locked pulses of the cr— polar- 
ized field, 

£:(t) =£:o(e-(*/*i)'-"-*+e-(*/*^)'-*('--^+-)*+'") . (1) 

For convenience, we use equal amplitudes. The energy of 
each state is denoted by e with an index for that state. 

The biexciton binding energy is A_| = 6+ + e_ — 

The dynamics of the four exciton levels is shown in 
Fig. |[ visualized as pseudospin (Bloch vector) for the 
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a— exciton in the absence or presence of the a+ exci- 
ton (see insert (a)). In the absence of cr+, the tr— ex- 
citon pseudospin evolves in an effective magnetic field 
Bx(i) = Bi?(i,fi)-|-Bij(t,i2), composed of a component 
Bi?- fixed in the a:— direction due to the resonant pulse 
with the |0} ^ |— ) transition and another B^, due to 
the off-rcsonant pulse with the 1+) ^ | H — ) transition, 
which rotates clockwise in the a;?/— plane at the beat fre- 
quency A+_. In the presence of the (t+, the pseudospin 
is driven by an effective Bxx(i) = Bfl,(t,ti) -I- Bi?(t, t2), 
in a rotating frame where the field component resonant 
with the biexciton resonance is fixed and the exciton field 
rotates anti-clockwise. The central idea is to optimize 
the control parameters (ti, 12) in such a way that, at 
the end of the shaped pulse, Bxx produces a tt rotation 
for the biexciton resonance, and Bx brings the exciton 
pseudospin back to the original state. We have studied 
numerically the dynamics of the four levels in the con- 
trol manifold generated by the parameters {ti, 12). Even 
using Rabi energies larger than the biexcitonic binding 
energy, it is possible to find a region in the phase space 
(^1,^2) fulfilling the required conditions. This is illus- 
trated in Fig. H(b) , where the evolution of the pseudospin 
in the absence and presence of the -I- exciton using a pulse 
shaped as in Eq. (|l|), with ti — 0.56 ps and ^2 = 1-05 ps, is 
represented on the Bloch sphere: at the end of the shaped 
pulse we have a tt rotation for |+) ^ | H — ), while |0} re- 
mains unaffected. As inferred from the calculation of the 
fidelity, this shaping works for an arbitrary initial state, 
and the pulse acts as a general rotor. A simple intuitive 
picture is given by the frequency spectrum of the pulse 
in Eq.(|I]), shown in Fig. |^(c). The two phase- locked com- 
ponents have opposite phases and the chosen ti and t2 
control the interference in the electric field giving a zero 
amplitude at the exciton resonance and a finite value at 
the biexciton resonance. The ultrafast manipulation of 
the qubits is carried out using light interference. For the 
7r/2 single qubit gate we used ti — t2 — 0.15 ps and in- 
phase components. The full theory for the application of 
the shaping to more general operations and limitations 
will be provided elsewhere [26| . 

A cogent illustration of the time limitation of deco- 
herence and the designs to beat the limit is given by a 
numerical simulation of a complete quantum computing 
algorithm. We choose the two-qubit DJ problem ||2l| ] 
both for its simplicity and for its fundamental nature. 
It has been implemented by nuclear magnetic resonance 
(NMR) For a single qubit variable x, there are four 
possible single-qubit functions: fi{x) — and f2{x) = 1 
belong to the class of constant functions and fz{x) = x 
and /^{x) = I — X the class of balanced functions. The 
DJ problem is to determine whether a given function is 
constant or balanced. Classically, the only way to find 
out if the function is constant or balanced is to calcu- 
late the function for both the and 1 input values and 
then compare the results. The quantum algorithm de- 



signs a unitary transformation which parallel-processes 
all possible values of x. We use the presence or absence 
of (T± excitons to represent the two qubits. The com- 
plete solution consists of three steps: (i) Encoding trans- 
forms the initial ground state dot into an input state: 
|m) = + -) + 1+) - |-) + |0)) by the single-qubit 
CT+ — 7r/2 rotation and ct— 7r/2 rotation, (ii) The uni- 
tary transformation Uf.\y,x) — Rx{TTy~^^'^^^\y,x), as- 
sociated to the given fj is applied. This is a conditional 
rotation acting on qubit x {a— exciton), conditioned on 
the effect of fj on the qubit y {a+ exciton). (iii) De- 
coding by 7r/2 rotations of both circular polarizations, 
ending with the ground state if the function is a con- 
stant and the |-|-) state if the function is balanced. The 
final measurement on the a+ exciton can be made with 
the cross-polarization probe or with the emitted photon 
polarization. 




FIG. 2. Dynamics of spectral shaping, (a) Energy schemat- 
ics for the transitions generated by a a— pulse, (b) Evolu- 
tion of the exciton sx and biexciton sxx pseudospins on the 
Bloch sphere under a shaped pulse, (c) Fourier transform of 
the shaped pulse and its components. 

The evolution of the qubits through the DJ algorithm 
is shown in Fig. ^ as the real parts of the coefficients of the 
corresponding states in two rows for a constant (/2) and a 
balanced (/a) function. The first column depicts the ideal 
case of sharp resonance functions without dephasing as 
described above. Even with the maximum packing of 
the pulses (using overlap at 10% of the peak electric field 
without substantial deterioration of the results) , the first 
column show the time to complete the quantum compu- 
tation above the adopted decoherence time of the order 
of 40 ps. Simple reduction of the pulse time to within 
dephasing time while keeping the pulse area leads to a 
failure to discriminate the solution, as shown in column 
2. Shrinking the temporal width of the pulses causes a 
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loss of the spectral selectivity needed to distinguish be- 
tween the excitonic and biexciton resonances. In order 
to exhibit the effects of the unintended dynamics and of 
the pulse-shaping remedy, the first three column have not 
included decoherence. The effect of shaping on quantum 
computation is shown the third column ]3C| ] : the possibil- 
ity to distinguish between constant and balanced within 
the decoherence time is recovered. 
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FIG. 3. Evolution of the qubits during a quantum compu- 
tation. The real parts of the coefficients TZcjit) for the first 
four levels are plotted. See text for explanation. 

The robustness of our scheme is tested with the in- 
clusion of decoherence in column 4. The pulse shaping 
can in principle be applied to more than two excitons. 
There is still the question of whether the pulse shaping 
is "scalable" to a system of n qubits, as n increases. The 
advantage of a quantum computer over a classical one is 
the saving of effort from 2" steps to order of powers in 
n Q . We need to avoid using an algorithm which treats 
the n excitons locked together as a strongly interacting 
system. The method is to process in steps, each in- 
volving a few excitons. Then, the overhead of shaping is 
at most a few powers of n [^, i.e., it does not spoil the 
exponential saving. 

In conclusion, we resolve the basic physics issues raised 
by the optical control of two antiparallel-spin excitons in 
a quantum dot. The application of this idea to the con- 
trol of two spin-polarized excitons in a single dot is well 
within present experimental capabilities, and we hope 
that our explicit design will stimulate a quick realization 
of an elementary quantum control. Our solution of pulse- 
shaping is a special case of quantum control of interact- 
ing spins and may be applied to other physical systems, 
such as trapped ions. Future work would involve exci- 
tons in different dots. Several schemes already exist for 
scalable system where qubits stored in an array of dots 
can be swapped around in a distributed quantum com- 
puter . Cavity Quantum Electrodynamics effects can 
be used for the interconnection between dots, and have. 



in principle, the additional potential of controlling the 
spontaneous emission of the excitons. 
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FIG. 4. Slowly varying envelopes of the pulse sequences 
used in the quantum computation for the solution of 
the DJ problem in Fig.^ The red (black) line indi- 
cates a+ (o" — ) circular polarization. The two rows la- 
beled by LUi refers to the resonance frequency of the pulses: 
LUi = e+ = e-,u)2 = e+ — A-| — = e_ A+_. 
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TABLE I. Pulse sequences used in the computa- 
tion, (1) and (2) indicate the two components of 
a shaped pulse of the form £{t) = £i{t) + £'2{t) 
= £:o(e-(*/*i)'-'"i*-"^i + e-('/t2)'-.-2t+»02) ^ The area 6 is 
defined as Ti9 = J d£{t) dt. dSo is 0.2 meV in the long pulses 
case and 2 meV for short and shaped pulses. The pulse se- 
quence of Uf-^{Uf^) is identical to the Uf^iUf^) except for a 
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